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Abstract 

Entropy production by the dilaton decay is studied in the model where the 
dilaton acquires potential via gaugino condensation in the hidden gauge group. 
Its effect on the Affleck-Dine baryogenesis is investigated with and without 
non-renormalizable terms in the potential. It is shown that the baryon asym- 
metry produced by this mechanism with the higher-dimensional terms is di- 
luted by the dilaton decay and can be regulated to the observed value. 



I. INTRODUCTION 

Many gauge-singlet scalar fields arise in the effective four- dimensional supergravity which 
could be derived from string theories. Among them the dilaton S has a flat potential in all 
orders in perturbation theory [1]. Therefore some non-perturbative effects are expected to 
generate the potential whose minimum corresponds to the vacuum expectation value (VEV). 
The most promising mechanism of the dilaton stabilization and the supersymmetry breaking 
is the gaugino condensation in the hidden gauge sector [2-5]. 

In cosmological consideration, even if the dilaton acquires potential through such non- 
perturbative effects, there are some difficulties to relax the dilaton to the correct minimum, 
as pointed out by Brustein and Steinhardt [6]. Since the potential generated by multiple 
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gaugino condensations is very steep in the small field-value region, the dilaton would have 
large kinetic energy and overshoot the potential maximum to run away to infinity. 

As a possible way to overcome this problem, Barreiro et al. [7] pointed out that the 
dilaton could slowly roll down to its minimum with a little kinetic energy due to large 
Hubble friction if the background fluid dominated the cosmic energy density. As a result, 
the dilaton could be trapped and would oscillate around the minimum. 

When the dilaton decays, however, the remaining energy density is transformed to radia- 
tion. One may worry about a huge entropy production by the decay, because it could dilute 
initial baryon asymmetry [8] . Therefore in order to obtain the observed baryon asymmetry, 
as required by e.g. nucleosynthesis, it is necessary to produce a larger asymmetry than the 
observed one at the outset. 

The attractive mechanism to produce large baryon asymmetry in supersymmetric models 
was proposed by Affleck and Dine [9]. However, the baryon-to-entropy ratio produced by 
this mechanism, n^/s, is usually too large. So if we take into account the entropy production 
after the baryogenesis, we can expect that the additional entropy may dilute the excessive 
baryon asymmetry to the observed value, as pointed out in e.g. [10,11]. 

In this paper, we investigate whether the dilution by the dilaton decay can regulate the 
large baryon asymmetry produced by the Affleck-Dine mechanism to the observed value. The 
paper is organized as follows. In §11 and §111, we describe the potential and the dynamics of 
the dilaton. Then, in §IV we estimate the baryon asymmetry generated by the Affleck-Dine 
mechanism taking into account dilution by the dilaton decay. Section V is devoted to the 
conclusion. We take units with 8nG = 1. 



We consider the potential of the dilaton non-perturbatively induced by multiple gaugino 
condensates. In string models, the tree level Kahler potential is given by 



where S is the dilaton, T is the modulus and $ represents some chiral matter fields [12]. 

The effective superpotential of the dilaton [2] generated by the gaugino condensation is 
given by 



II. DILATON POTENTIAL 



K = — \n(S + S*) - 3 ln(T + T* - |$| 2 ), 



(1) 



—a a S 



(2) 
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Here, for the SU(N a ) gauge group and the chiral matter in M a (N a + N a ) "quark" represen- 
tations, a a and A a are given by 



Then the modulus T have a potential minimum due to the presence of the Dedekind function 
i](T) [5]. Hereafter we assume the stabilization of the modulus T at the potential minimum 
and concentrate on the evolution of the real part of the dilaton field. 

Since at least two condensates are required to form the potential minimum, we consider 
a model with two condensates. Then the indices of the gauge group are a = 1,2, and we 
take 10 ^ «i ^ «2- 

The potential V for scalar components in supergravity is given by 



V = e K 

where 



{K- L ) t jD i W{D j Wy - Z\W\ Z , (5) 



dW dK 

Kj = d 2 K/d¥d$*, the inverse (K' 1 )^ is defined by (K^YjK^ = 8\ and i = S,T, $. In 
the region, aReS* ^> 1, the potential Eq. (5) can be rewritten as 



V{S)-e K {K' 1 ) s s \d s W\ 2 

= (,S + 1 S*)|aiA 1 |V ai(5+5 * ) 



1 + ?^l e -{*2- ai )S 2 



(7) 



In this potential the imaginary part of S has a minimum at 



tafiU. = <^i±^, (8) 

where n is an integer. So we assume ImS* = ImS^n and concentrate on the behavior of KeS 
hereafter. Then we find that the potential minimum, ReS'min, is given by 

ReS min = ^— ln(^V (9) 



Q.2 — Oil \Ct\ Ai 

We assume ReS'min — 2 to reproduce a phenomenologically viable value of the gauge coupling 
constant of grand unified theory [5]. The negative vacuum energy at the minimum of the 
potential is assumed to be canceled by some mechanism such as a vacuum expectation value 
of three form field strength. Although we consider a model with two gaugino condensates, 
our following estimations would be almost unchanged in a single condensate model with non- 
perturbative Kahler corrections [13,14], because the evolution of the dilaton is determined 
only by the slope of the potential in the region S S min and its mass. 

On the other hand, the position of the local maximum of potential, ReS^ax, is given by 

ReS max = ReS min + In (^) . (10) 



«2 — Oil \Cti . 

For S <^ S min , the potential (7) can be approximated as 

V(S) ~ 2ReSV e~ a22ReS , (11) 
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where Vb — |«2A2| 2 . The potential (7) has the minimum at S m i n and the local maximum 
at S max (> S m [ n ). For S ^, S CT = S min — l/(a 2 — eei), the approximate expression for the 
potential (11) breaks down. Around the potential minimum S min , the potential (7) becomes 



V(S) ~ \a ± A 



-2ai5 n 



(a 2 - aO^ReS" - ReS min ) 2 . 



(12) 



However, one can see from the Kahler potential (1) that the variable ReS* does not have 
the canonical kinetic term. Therefore we introduce the canonically normalized variable, 0, 

as 



= —= In ReS. 
V2 



(13) 



III. STABILIZATION MECHANISM FOR THE DILATON 



Here, after reviewing the mechanism for dilaton stabilization proposed by Barreiro et 
al. [7], we estimate the relic energy density of the dilaton and the amount of the entropy 
density produced by its decay. We will consider the situation that the universe after inflation 
contains the dilaton, 0, and a fluid with the equation of state, p — (7 — l)p, where 7 is a 
constant. For example, 7 = 4/3 for radiation or 7 = 1 for non-relativistic matter. The latter 
includes oscillating inflaton field or/and the Affleck-Dine (AD) condensate 4>ad- 

In the spatially flat Robertson- Walker space-time, 



ds 2 = -dt 2 + a(t) 2 dx 2 , 
with the scale factor a(t), the Friedmann equations and the field equation for <f> read 

H 

2 " 

dV(4>) 



<p = -3H<p 



d<f) 



H 2 



1 • 



P 



r + v{<i>) 



(14) 

(15) 
(16) 
(17) 



where H = a/a is the Hubble parameter and a dot denotes time differentiation. We define 
the new variables, 



x = 







V6H' 



y = 



(18) 



and the number of e-folds iV = ln(a). 

Then, the equations of motion can be rewritten as 



x 

y' 

H' 



-3x - y - 



2v y+ r 



2x 2 + 7(1 - x 2 - y 2 ) 



, , ■ xy + -y [2x 2 + 7(1 - x 2 - y 2 ) 



3 -H 

2 



2x 2 + 7(l 



x 



y 2 ) 



(19) 

(20) 
(21) 
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where the prime denotes a derivative with respect to N. In terms of these variables, the 
Friedman equation becomes x 2 + y 2 + pj (3H 2 ) = 1. We see that x 2 and y 2 are respectively 
the ratios of the kinetic and potential energy densities of the dilaton to the total energy 
density. We consider the case of the universe dominated by the background fluid, so the 
inequalities, x 2 ,y 2 <C 1, hold. Then Eq. (21) can easily be solved and the solution is 



H = H e 



-37JV/2 



(22) 



Next we introduce another new variable, 



_ dReS 
x s = — x. 

dip 

Then Eqs. (19) and (20) can be respectively rewritten as 

. 3 [3 fdReS\ 2 2 

*s = -^s + -jx s + 2a 2] j-^—) y, 

/3 3 
2 x »y + 2 1V ' 



(23) 



(24) 
(25) 



where we have used the relation —2a{V = dV/dReS. Now we examine the stationary points 
in the above equations. From y' — 0, we find 



x s 




(26) 



On the other hand, from x' s = we see 



2 3(2-7)7 f dReS \ 
V 2{2a 2 ) 2 \ d<P J 



(27) 



Except for the factor (dReS/d<p)~ 2 , Eqs. (26) and (27) represent the scaling solution for 
the scalar field with exponential potential [15]. In spite of the presence of the factor 
(dReS '/ d(j))~ 2 , we can verify that the deviation from the scaling solution is small enough, as 
already shown in [7]. Thus we can write the solution as 



37 1 
ReS = — ^-iV + — In 



2a 2 



a 2 



(2-7)7 V3#o 



2a, 



+ e(N), 



(28) 



where 



e(N) = — In N + — In 
a 2 1 2a 2 a 2 



4K, 



(2-7)7 V3#o 



2a, 



(29) 



denotes the deviation from the scaling solution; we find x' s = e(N)" <^ 1. 

Figure 1 depicts time evolution of ReS as a function of N for various initial values of the 
dilaton field amplitude. As is seen there if the dilaton relaxes to the scaling solution before 
reaching S m [ n , its energy is small enough to prevent overshooting. Attractor behavior in a 
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different situation has been studied in [16]. The Hubble parameter at S — S m i n is estimated 

as 



2a 2 2V Smin x , . 

^nun- 3 y (2 _ 7)7 e -tnun, ^Uj 

from Eqs. (22) and (28). On the other hand, as is seen from Eq. (7), the mass of the dilaton 
in vacuum is given by 



~ 2{a 2 - a^^Voe-^^. (31) 

We therefore find ~ H min , so the dilaton begins to oscillate immediately when it ap- 
proaches the potential minimum. Since the mass of the gravitino is given by m 3 / 2 — 
A 2 e~ Q2,Smin , the mass of the dilaton is 

~ alm 3/2 

When the dilaton S approaches the critical point S cr , the single exponential approximation 
(11) breaks down and the scaling behavior terminates. Then the energy density of the 
dilaton is estimated as 



V 



3H 2 (x 2 + y 2 ) 




Indeed, the numerical calculation gives a close value ~ 10 _4 p at the beginning of the 
oscillation. After that, the dilaton begins to oscillate and the energy density decreases as 
a(t)~ 3 until it decays. 

In Fig. 2, we present the energy density of the dilaton in the universe filled by the 7 = 4/3 
background fluid at the beginning of the oscillation regime for various model parameters, 
«i and «2- The values along the contour lines represent the energy density p^ in the unit 
of 10 _4 p. The case with 7 = 1 is depicted in Fig. 3, where we find smaller energy density of 
the dilaton by a factor of ~ 3. These figures are drawn in the two-parameter space, though 
the potential contains four parameters as seen from Eq. (7). The other two have been fixed 
by setting m 3 / 2 = 1 TeV and S inin = 2 [5]. 

The decay of the dilaton produces huge entropy. If we assume that the background fluid 
is radiation with 7 = 4/3, at the time t = t m - m its temperature is T ~ 10 11 GeV and the 
entropy density is 

s = ^ 9 ,T 3 , (34) 
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where g* ~ 10 2 is the effective number of relativistic degrees of freedom. By using Eq. (33), 
we find that the entropy density increases by the factor 

~~ To p , " 10 U^GevJ I T D J I 10"* ' (35) 
when the dilaton decays at 

// = r^m^lO--Gev( I ^ v ) 3 , (36) 

where 

r D ^/^10->GeVy^) 3/ \ (37) 
is the reheating temperature after decay of the dilaton. 

IV. AFFLECK-DINE BARYOGENESIS 

The Affleck-Dine mechanism is an efficient mechanism of baryogenesis in supersymmetric 
models [9,17]. In fact, it is too efficient and the produced baryon asymmetry, n^/s, is in 
general too large. However, additional entropy release by the dilaton decay may significantly 
dilute the baryon asymmetry [10,11] and we examine this possibility here. 

It is known that the Q-ball formation occurs for many Affleck-Dine flat directions [17-19]. 
Whether Q-balls form or not depends on the shape of radiative correction to the flat direction 
[19,20]. In this paper we estimate the baryon asymmetry providing that the Affleck-Dine 
field does not lead the Q-ball formation. For instance, one example is a flat direction with 
large mixtures of stops in the case of light gaugino masses, another is i/ u L-direction [20]. 

A. Original Affleck-Dine mechanism 

First, we investigate the originally proposed Affleck-Dine mechanism with a flat potential 
up to (f>AD ~ 1 [9]. We consider the situation that there are the dilaton and the Affleck- Dine 
condensate in radiation dominated universe. As mentioned above, at the moment H = m^, 
the dilaton begins to oscillate with the initial energy density p^ ~ 10~ 4 p 7 , where p 7 is the 
energy density of the background radiation. Then, on the other hand, the AD condensate 
is expected to take a large expectation value, (ft ad ~ 1, above which its potential blows up 
exponentially. 

The amplitude of the AD field and its energy density remains constant while the Hubble 
parameter is larger than itlad, where tuad is the mass of the AD condensate. One should 
keep in mind, however, that though the energy density, Pad, remains constant the baryon 
number density decreases as 1/a 3 if baryonic charge is conserved. So if baryon charge is 
accumulated in "kinetic" motion of the phase of (ft ad it would decrease as 1/a 3 . If, on the 
other hand, the AD-field is frozen at the slope of not spherically symmetric potential then 
baryonic charge of the AD-field is not conserved and after H < itlad both radial and angular 
degrees of freedom would be "defrosted" and baryonic charge may be large. 
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As we noted above, when H ~ m 3 / 2 (— toad), the AD field begins to oscillate. Its energy 
density at that moment becomes comparable to that of the radiation, while the energy 
density of the dilaton is estimated as 



PAh-- 



1/2 



P± 
P 7 



Pi\h= 



m 3/2 



Py\H= 



™3/2 



(38) 



for the initial energy density = 10 _4 p 7 and ~ 10 2 toad- After that the universe 
becomes dominated by the oscillating AD condensate and enters into approximately matter 
dominated regime (see below). 

The energy density of the condensate and its baryon number density are given respec- 
tively by the expressions: 



Pad = m 2 AD (f) 2 AD , n b = Km AD 4> 2 ALn 



(39) 



where k = rib) 'n ad < 1 is a numerical coefficient and had is the number density of the AD 
field. 

The rate of evaporation of the condensate, given by the decay width of the AD field into 
fermions, T AD = Cm ad with C — 0.1 — 0.01, is quite large. When the Hubble parameter 
becomes smaller than Tad, thermal equilibrium would be established rather soon. However, 
the condensate would evaporate very slowly and disappear much later [21]. The low evap- 
oration rate is related to a large baryonic charge and relatively small energy density of the 
condensate. Below we will find the temperature and the moment of the condensate evapo- 
ration repeating the arguments of ref. [21]. Let us assume that the condensate evaporated 
immediately when H = Tad producing plasma of relativistic particles with temperature 
Tad and chemical potential Had- The temperature can be estimated as Tad — Pad an d 
since Tad m AD the chemical potential is given by 



Pad 



T 2 

1 AD 



k4>ad > m A D, 



(40) 



if k is not very small. On the other hand, chemical potential of bosons cannot exceed their 
mass. It means that instantaneous evaporation of the condensate is impossible. The process 
of evaporation proceeds rather slowly with an almost constant temperature of the created 
relativistic plasma. During the process of evaporation the energy density of the latter was 
small in comparison with the energy density of the condensate, except for the final stage 
when the condensate disappeared. 

The cosmological baryon number density and energy densities are given by the equilib- 
rium expressions: 



rib,tot 
T 3 

Ptot 
^4 



2N f N c B q 
6tt 2 



fe 3 + - 2 e 



1 



7r 2 g« 7 2N f (N c + l)n 
30 8 15 



?/ ' 2n 2 Jo 

2 



1 



1 



exp(e-£)-l exp(e + f)-l 

2 , r- / \ 4' 



+ B c (41) 



7T 



7T 



1 f°° 

2^ Jo dVr]2 



1 



exp(e-£)-l exp(e + 0~l 



(42) 
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where rj = p/T is dimensionless momenta, e = \Jv 2 + m AD/T 2 , £ = Pad/T and £ g = /x ? /T 
are dimensionless chemical potential of the AD field and quarks, respectively. Nf = 6 
and N c = 3 are the numbers of flavors and colors and factor 2 came from counting spin 
states, B q = 1/3 is the baryonic charge of quarks while the baryonic charge of the AD field 
is assumed to be 1, and B c and p c are baryon number density and energy density of the 
condensate normalized to T 3 and T 4 , respectively The first term in p tot includes energy 
density of light particles with zero charge asymmetry and g* is the number of their species. 
The second term includes the contribution from leptons with the same chemical potential 
as quarks - it is given by (N c + 1). 

For definiteness, let us assume that the AD field decays into the channel 4> AD — > 3q + / 
and taking into account that the sum of baryonic and leptonic charges is conserved 1 , so that 
B - L = 0, we find 

PAD ,. Q , 
Pq = Pl = —■ (43) 

Before complete evaporation of the condensate, the chemical potential of AD-field remains 
constant and equal to its maximum allowed value trad- Thus the only unknowns in these 
expressions are the temperature and the amplitude of the field in the condensate. According 
to Eq.(39), from Eqs. (41) and (42) we obtain 

^ = k — — — . (44) 

Pc m AD 

The same relation was true for the initial values of n^tot/^ 3 and p to t/^ 4 - Assuming that 
the ratio n^tot/ Ptot remains the same during almost all process of evaporation, though it is 
not exactly so, we can exclude B c ,nb ttot , p c and p tot from expressions (41) and (42) and find 
one equation that permits to calculate the plasma temperature in presence of evaporating 
condensate as a function of the baryonic charge fraction in the initial condensate, k. We 
find 

m AD /T ~ 20, for « = 1, (45) 
m AD /T ~ 2, for k = 0.1. (46) 

Exact solution of the problem demands much more complicated study of the evolution of 
energy density according to the equation p = —3H(p + P), while the evolution of baryonic 
charge density is determined by the conservation of baryonic charge which is assumed to 
be true at the stage under consideration and thus n^tot oc a~ 3 . The temperature of plasma 
found in this way would not be much different from the approximate expressions presented 
above. 



One may wonder if this assumption is inappropriate because the baryon asymmetry created in 
this channel would be washed out by anomalous electroweak processes [22]. As will be seen later, 
however, we can avoid this difficulty because in most cases of our interest the AD condensate 
evaporates at a lower temperature when these anomalous processes are no longer effective. 
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Using the above-calculated plasma temperature (46) we find that at the moment of 
condensate evaporation (when pad — m AD) the cosmological energy and baryon number 
densities of the created relativistic plasma are given by 

p p ~ 1000T 4 , (47) 

n b ~ 50T 3 , (48) 

for k — 1, and 

P P ^ 70T 4 , (49) 

n b ~ 1.75T 3 , (50) 

for k = 0.1. 

We see that a large baryon asymmetry prevents from fast condensate evaporation, though 

the interaction rate could be much larger than the expansion rate. From Eq. (48) or (50) 
and the baryon number conservation 

Tib — KITlAD ( —JW ) {4>AD\H=m AD ) 2 , (51) 



we find 



a ev \ 3 m A D 



) = K ((f>AD\H=rn AD ) 2 (52) 




2 



<t>AD\H=m AD \ ^ 1Q 3 3) for K=1> (53) 



m AD 



2 



046 {^Lff=m^j ^ 3xl0 3 0) for K = Q-1 (54) 
m AD J 

where a^z) and a ev are the value of the scale factor at the moment H = trad and that at the 
evaporation of the AD field, respectively. Then at the evaporation the Hubble parameter 
and the baryon-to-entropy ratio are respectively given by 

H - (o /3)^ - { 2 X 10_14 GeV for K = f55) 
n cv - [ Pp /6) - 5 x 10 _i3 GeV for K = QA W 



nb 
s 

rib 
s 



ev 



1, for k = 1, (56) 
-0.04, for « = 0.1, (57) 



from Eqs. (47), (48), (49) and (50). 

Now we have to calculate the ratio of the baryon asymmetry to the entropy of the plasma 
after thermalization of the products of dilaton decay. Initially, at the moment of evaporation 
of AD-condensate the energy density of the dilaton is roughly 10~ 3 with respect to the energy 
density of plasma. The latter is dominated by chemical potential \x = vtiad > T. When the 
universe expanded by the factor Pad/P4>\ cv = o C q/a C v — 10 3 the dilaton starts to dominate 
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and the relativistic expansion regime turns into matter dominated one at a = a eq . To the 
moment of the dilaton decay the energy density of the dilaton becomes larger than the 
energy density of the plasma formed by the evaporation of the AD-condensate by the factor 
Od/ocq = (H cq /H d ) 2 / 3 , where a d and H d are the scale factor and the Hubble parameter at 
the time of the dilaton decay. Keeping in mind that H eq = (a ccl /a cv ) 2 H cv = 10~ 6 H ev , we 
obtain the dilution factor by the dilaton decay 



3/4 



,Pad 



( H cv al 



1/2 



Pad 



H d 



1/2 



(58) 



Thus for k — \ the dilution factor is only 14, while for k — 0.1 it is 70. 

Finally we find that the baryon asymmetry after the dilaton decay is given by 



s 



s 



— ~ 0.1 

cv A 



0.001. 



(59) 



Thus in this model the dilution of originally produced asymmetry from the decay of AD 
field is not sufficient. 



B. Affleck-Dine mechanism with non-renormalizable potential 

As we have seen, additional entropy production due to the dilaton decay is too small 
to dilute the baryon asymmetry generated in the original Affleck-Dine scenario. Therefore 
it is necessary to suppress the generated baryon asymmetry. The presence of the non- 
renormalizable terms can reduce the expectation value of the AD field during inflation. As a 
result, the magnitude of the baryon asymmetry can be suppressed. Hence we introduce the 
following non-renormalizable term in the superpotential to lift the Affleck-Dine flat direction 
as a cure to regulate the baryon asymmetry [23,24]. 

W = * « <Fad, (6°) 

where M is some large mass scale. 

The potential for the AD field in the inflaton-dominated stage reads 

V( M = + (^§^ + H.C.) + |A|'^t!, (61) 

where c\ and C2 are constants of order unity. The first and the second terms are soft terms 
which arise from the supersymmetry breaking effect due to the vacuum energy of the inflaton. 
The minimum of the potential is reached at 

1^1 = (™_) . (62) 

During inflation, the AD field takes the expectation value \4>ad\ — (H in{ M n ~ 3 / \y^ n ~ 2 \ 
where H in{ is the Hubble parameter during inflation. After inflation, it also traces the 
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instantaneous minimum, Eq. (62), until the potential is modified and the field becomes 
unstable there. The AD field starts oscillation when its effective mass becomes larger than 
the Hubble parameter. There are three possible contributions to trigger this oscillation: the 
low energy supersymmetry breaking terms, the thermal mass term from a one-loop effect 
[25] , and the thermal effect at the two-loop order [26] . 

First we consider the case that the low energy supersymmetry breaking terms are most 
important and that these thermal effects are negligible. When H ~ "023/2, the low energy 
supersymmetry breaking terms appear and the potential for the AD field becomes 

V(M = m * AD M + + H.C.) + |Ar^£!, (63) 

where A is a constant of order unity, and the ratio of the energy density of the AD field pad 
to that of the inflaton pj is given by 

PAR „ ( rn 3/2 M \ (g4) 

up to numerical coefficients depending on ci, C2 and A. For example, the typical value for 
n = 4 is pad/ Pi ^ 1(T 16 (M/A). 

Let us first assume that the inflaton decayed at H — rriAD(— ^3/2) with the decay rate 
Tj = rriAD and after that the universe was dominated by radiation. Note that the corre- 
sponding reheating temperature is T R ~ ^/T] ~ -Jtoad — 10 10 GeV. Then the evaporation 
of the AD condensate into relativistic plasma would be different from the evaporation into 
cold plasma considered in the previous subsection. Due to interaction with plasma the prod- 
ucts of the evaporation acquire much larger temperature than in the case of the evaporation 
into vacuum. Since the energy density of the condensate is negligible in comparison with 
the total energy density of the plasma, the temperature of the latter drops in the usual way, 
T oc 1/a, in contrast to the previously considered case when T = const. Since the temper- 
ature of the plasma is high, T ^> rriAD, the baryon number density is rib = B C T 3 + CsT 2 p 
where Cb ~ 1 is a constant coefficient, p < ttiad is the value of the chemical potential, and 
we have neglected terms of the order of p 3 . Since n b oc a -3 , the ratio of a cv to a ad is 



n b \ H -- 



■m A D 



O-AD TTIAdTb 



(65) 
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( ^)(~y ( M )fo _ 4 , 



(compare to Eq. (52)). Here we took for the initial value of the baryonic charge density 
rib\H=m AD = kijiad^ad w hh k ~ 1. For n — 4, a cv /aAD becomes of order 10, and we find 
that the condensate would evaporate soon. 

To be more precise, however, we must take into account that the interaction rate 
of the condensate is Tad = (0.1 — 0.01)mAD and the evaporation cannot start before 
H = Tad- At that moment the plasma temperature would be smaller by the factor 
(m ad/Tad) 2 = 10 2 — 10 4 and the baryon number density of the condensate would be smaller 
by (t^ad/Tad) 6 - Correspondingly the red-shift of the end of evaporation should be shifted 
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by a factor (ttiad/^ad) 2 with respect to the beginning of evaporation and it means that it 
would remain the same with respect to the initial moment H = itlad- 

As we have already noted, for n = 4 the condensate decays quickly and the baryon 
number density produced in the decay is diluted by the plasma created by the inflaton 
decay as 

nj, T^P^o = 10 _ 9 t T R \ ZlTeVN ,M\ 

s m AD p, \WGeV)\m AD ) \ X) K ' 

where Tr is the reheating temperature of the inflaton and we used the estimate of Eq. (64). 
The result does not depend upon the moment of the decay of AD-condensate, since its 
energy density remains sub-dominant. If an additional dilution by the dilaton and the early 
oscillation by a thermal effect [25] are operative, the baryon asymmetry become even smaller 
than the observed one and the n = 4 model cannot explain the observed baryon asymmetry. 
Hence we must consider the flat direction with n > 4. 

Hereafter, we study the AD fields with n > 4 including the thermal effect. For the AD 
fields with n > 4, the relevant thermal effect comes from the running of the gauge coupling 
constant [26] rather than the thermal plasma effect [25]. The potential for the AD field in 
the inflaton-dominated stage reads 



2n-2 

(68) 



V^ad) = (- Cl # 2 + m^)|0 AD | 2 + aT 4 m 

| ( c 2 \H(j) n AD | Am 3/2 (j)AD lEC \ | \\\^ AD 




nM n ~ 3 nM n ~ 3 j M 2n - Q 

where the second term is the thermal effect at two loop order which is pointed out in [26] 
and a denotes the gauge coupling. 

If the effective mass of the AD field becomes comparable to the Hubble parameter when 
it is larger than the low energy supersymmetry breaking scale, 

~|2 



a-^—^H 2 (>m AD ), (69) 



>AD\ 



then AD field undergoes the early oscillation by the thermal effect. During the oscillating 
inflaton dominated stage (t < t r h), the temperature of the plasma behaves as 



R 



a{t lh ) \ ^ T i,2 H w (70) 



a(t) J 

From Eqs. (62) and (70), the effective mass term of the AD field is rewritten as 



T 4 / ^ \2/("-2) 



a- 



>AD\ 



12 



-^(^1 H^l^\ (71) 



By comparing with Eq. (69), when the AD field begins to oscillate at t = t os , the Hubble 
parameter is given by 

/ \ (n-2)/n ( \ \ 2 ^ n , , 

K«) ' (72) 
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From now on we concentrate on the case n = 6, for which 




H os ~ ITeV I -7—^ ) TT^ry TFf ■ ( 73 ) 



Thus we find that the AD field begins to oscillate at if ^ ^ad due to the thermal term if 
> 10 8 GeV for M = 1, and if > 10 6 GeV for M = 10~ 2 , respectively. During the 
early oscillation driven by the thermal term, the amplitude of the AD field decreases as 

/ \ 9/4 

|0AD(t)| = |0ADL^J , (74) 

where a os denotes the scale factor at the beginning of the oscillation. The analytic deriva- 
tion of Eq. (74) is shown in Appendix and we have confirmed this result by the numerical 
calculation. It also agrees with the analysis in [27]. Using Eqs. (62), (70) and (74), we 
obtain the ratio of the amplitude of the AD field to temperature of plasma as 

*ad\ ( aos \ 15/8 

(75) 



T \t*\) {a(t)) 



We consider the evaporation rate of the AD field. The condition that the particles coupled 
to the AD field are light enough to exist as much as radiation reads < T, where h is 

the corresponding coupling constant. Following [25] let us adopt the scattering rate of the 
AD field T ~ h 4 T in this situation as the rate of its evaporation. Estimating }i\4>ad\/T and 
T at the reheating time, we find 



h\4>. 



AD 



10 -2\ 5 /6 




T 
and 

T h 4 /l0 10 GeV\ / h 




(76) 



H T R \ T R / V 10 



-2 



(77) 



Eq. (76) shows that the particles coupled to the AD field are thermally excited and populated 
well before the reheating time for any reheating temperature, T R < 10 16 GeV, while we find 
that the AD condensate can evaporate around the typical reheating time from Eq. (77). 
The baryon number density for the AD field (Pad is given as 

n b = -iq((f)* AD <j) AD - (f) AD (f) AD ), (78) 

where q is a baryonic charge for the AD field. 

The baryon number density at H = H os is estimated as 



_ 4gm 3 / 2 6 , 



Aq5m 3/2 (H os M*\ 1/2 



3A V A 



where 5 is a effective relative CP phase. The baryon-to-entropy ratio at the reheating time 
is estimated as 



14 



s 



3Tr n b 



Since the dilution factor by the dilaton decay is given by 



qSm 3/2 T R { H os M* \ 1/2 

— H[ s {~^) ' (80) 



' m 1 n 



10 4 T D V 10- 

from Eqs. (80) and (81), we obtain the final baryon asymmetry, 



(81) 



rh _ nb 
s s 



1 

A 



in _ 8 r ^s/2\ 3 / 2 /M\ 3 / 2 / T g \ ZlTeV\ 




(82) 



for if os > 7713/2. This result can easily meet the observation if we take, for example, H os ~ 
10 2 m 3 / 2 and other factors to be of order of unity. In this case, from Eq. (77), we find that for 
the present case the AD condensate can evaporate before the inflaton decay is completed. 

On the other hand, in the case the reheating temperature is so low that the thermal 
effect does not lead the early oscillation, the reduction of the cut-off scale M can lead to the 
reasonable baryon asymmetry: 



H± _ in -n. 
s 



This expression implies that the cut-off scale M should be around the GUT scale 10~ 2 or 
10 16 GeV for inflation models with a low reheating temperature, Tr < 10 6 GeV. Furthermore 
we find that the final baryon asymmetry is independent of the reheating temperature of 
inflation within the range. 

In the present model, the supersymmetry breaking is caused by the F-term of the dilaton. 
Therefore, when the dilaton decays, it can decay into gravitinos through their mass term 
and this process could lead to overproduction of the gravitinos. The constraint derived in 
[28] to avoid the overproduction is > 100 TeV. The mass of the dilaton, (32), in the 
model considered here is in the allowed region. 



V. CONCLUSION 

In this paper, we have studied the Affleck-Dine baryogenesis in the framework of the 
string cosmology. In string models, the dilaton is ubiquitous and does not have any potential 
perturbatively. We adopted the non-perturbatively induced potential of the dilaton via the 
gaugino condensation in the hidden gauge sector. Then we set phenomenologically desired 
values for the gravitino mass and the VEV of the dilaton. 

The attractive mechanism to stabilize the dilaton at the desired minimum was proposed 
by Barreiro et al. [7]. They did not estimate the energy density of the oscillating dilaton. 
It is estimated in the presented paper where we have found ~ 10 _4 p at H = m^. This 
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energy transforms into the radiation after the decay of the dilaton before nucleosynthesis 
because the mass ~ 10 2 TeV is sufficiently high. 

We have discussed cosmological baryogenesis in this model. In the above-mentioned cos- 
mological history with the entropy production, the Affleck-Dine baryogenesis might be the 
only workable mechanism for baryogenesis. We have investigated the Affleck-Dine baryoge- 
nesis with and without non-renormalizable terms. We have shown that while the original 
Affleck-Dine scenario produces too much baryon asymmetry even if there is the dilution by 
the dilaton decay, the model with n = 6 non-renormalizable terms can lead to the appropri- 
ate baryon asymmetry. 
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APPENDIX 

In this appendix, we derive Eq. (74). Although the similar discussion can be found in 
the literature [29], we derive it for completeness. 

We consider the evolution of the AD field after the beginning of oscillation induced by 
the thermal effect at the two-loop level. Then the AD field obeys the equation of motion 

4>ad + 3H<j) AD + a— = 0. (Al) 

By decomposing <Pad into 

0ad = \4> AD \e id = <&e* e , (A2) 
Eq. (Al) is reduced to the following equations 

<t> + 3#<i> - 6 2 § + a— = 0, (A3) 

(a 3 9<5> 2 } = 0. (A4) 

The second equation (A4) is interpreted as the conservation of the angular momentum which 
corresponds to the baryon number density and rewritten as 



6& = 6& 



a 



3 

os \ ^2 



= m*t\^r\ , (A5) 



a(t) J \a(t) 



where <3>o represents the initial amplitude of the AD field and m means the initial angular 
velocity of the order of m 3 / 2 for n — 6. By eliminating 6 in Eqs. (A3) and (A4), we obtain 



2 f a„\° fi„\' _ T 4 



• +aH *- mJ (^J ItJ * +a T = °- (A6) 
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Multiplied by $ and using T 4 oc cr 3/2 , Eq. (A6) yields 



d 
~dt 



A2 , 2 / u os 

<P + m 



6 ^ $ 2 



a(t) / $ 2 T 2 



= -QH 



$2 + m 2^«os^ ^ 



S CP 2 S 

|tfaT 4 ln- + ^aT 4 . 



a(t)/ $ 2 

On the other hand, multiplying Eq. (A6) by $, we obtain 
1 



(A7) 



6 $4 



a(t)J $ 2 



^ + a r 4 = 0. 



(A8) 



By taking the time average over the time scale of the cosmic expansion, we obtain the cosmic 
virial theorem 



a(t)J $ 2 



(A9) 



where (...) denotes the time average. Moreover, since the second term which represents the 
centrifugal force becomes efficient around only $«0, Eq. (A9) could be rewritten as 



(cb 2 ) = ( a T A ). 
From Eqs. (A7) and (A10), we obtain 

d (, 1 $2 \ 15 rr 



Since we know T oc a 3 / 8 , hence, we find 



(A10) 



(All) 



(A12) 



for $ > T. 
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FIGURES 




FIG. 1. Evolution of Re S as a function of N for various initial conditions with SVmn = 
in case of 7 = 4/3. We set the initial values of the velocity and the Hubble expansion rate 
dReS/dt\o = 0, and Hq = 1, respectively. 
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FIG. 2. Energy density of the dilaton at the beginning of the oscillations for various model 
parameters. The number associated with each contour line represents the value of p^ normalized 
by 10 p. Here we adopt 7 = 4/3. 
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FIG. 3. Same as Fig. 2 except for 7 = 1. 
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